PERSISTENCE OF GAPS IN THE SPECTRUM OF CERTAIN 
ALMOST PERIODIC OPERATORS 
Norbert Riedel 

Abstract: It is shown that for any irrational rotation number and any admissible 
gap labelling number the almost Mathieu operator (also known as Harper's 
operator) has a gap in its spectrum with that labelling number. This answers 
the strong version of the so-called "Ten Martini Problem" . When specialized to 
the particular case where the coupling constant is equal to one, it follows that 
the "Hofstadter butterfly" has for any quantum Hall conductance the exact 
number of components prescribed by the recursive scheme to build this fractal 
structure. 



Introduction 

The present work is concerned with the spectral properties of the simplest kind 
of discrete Schroedinger type operators with an almost periodic potential. These 
operators form a self-dual class with respect to the Fourier transform, with ex- 
actly one operator being invariant. More specifically, a precise one-to-one rela- 
tionship between gaps and admissible gap labelling numbers will be established. 
While relegating the technical formulation of the problem and the outline of its 
proof to the first section, we are going to dwell for the rest of the introduction 
on the significance of this result in the special case of the self-dual operator for 
physics and mathematics. 

Recent progress that has occurred in solid state physics through the development 
of new experiments which led to improved measurements, has made it possible 
to obtain improved physical evidence for the presence of the butterfly fractal 
spectrum. In [GSUKNKS] portions of the "Hofstadter butterfly" have been 
observed in lateral superlattices patterned on GaAs/AlGaAs heterostructure 
by exploiting the quantum Hall conductance as a diagnostic tool. One way to 
view the significance of the present paper is the recognition, that this diagnostic 
approach to detecting the fractal stucture experimentally is on solid theoretical 
ground, at least as long as the "Hofstadter butterfly" is accepted as a paradigm 
for the quantum Hall effect: Knowing the specific value of the quantum Hall 
conductance, it is possible to allocate the corresponding gap components in the 
fractal structure in a prescribed way, without running the risk of missing one. 
This follows from the considerations in [OA], where it was shown that, assuming 
that the strong version of the "Ten Martini problem" holds, the components 
associated with the same quantum Hall conductance can be counted by means 
of a specific combinatorial formula. In another experiment cold neutral atoms 
in an optical lattice were used to exhibit salient features of the butterfly fractal 
[JZ] . The employment of lasers allows the simulation of a magnetic flux through 
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the lattice. The success of this approach for specific rotation numbers depends 
on the visibility of the particle density, which in turn depends on its periodicity. 
As shown by a specific experiment in [JZ] (FIG 4b), the visibility decreases for 
the irrational rotation number a — ^ . Again, as the present work will show, it 
is reassuring to know, that the butterfly fractal does not allow for unexpected 
transitions at irrational rotation numbers to occur, which could not be picked 
up through a suitably designed experiment. 

We turn now to the mathematical significance of the present work, which at 
this point is somewhat speculative in nature. The set of quantum Hall numbers 
is a cyclic subgroup of the additive group of real numbers. If g is the posiitve 
generator of this group, then each quantum Hall number q can be written as q = 
kg, for a suitable integer k. As was shown in [OA] under the assumption that the 
strong version of the "Ten Martini conjecture" holds, the number of components 
in the butterfly fractal with a common positive quantum Hall number q is equal 
to 

n 

$(2fc), where = ^ <f(j), 

and (p is Euler's totient function. On the other hand, an observation by J. Franel 
from the 1920's asserts that the Riemann hypothesis is true if and only if the 
relation 

(n) 

holds for any e > 0. Here is the jth Farey fraction of order n. For a detailed 
exposition of this subject see Edmund Landau's lectures on number theory [L], 
Band 2, Kapitel 13. Notice that the Farey fractions, in the order as they appear 
in this asymptotic formula, occupy a natural position in the butterfly fractal. 
So it appears that the butterfly fractal holds information about the Riemann 
hypothesis, the exact nature of which remains yet to be determined. Inspired 
by the color coding of the butterfly fractal according to the Hall conductance, 
which was also introduced in [OA] , and which has become very popular in recent 
years, it is tempting to take a cue from Marc Kac, who is not only alleged to 
have offered ten martinis as a reward for the solution of the eponymous problem 
which is the subject of this work, but who also famously asked "Can one hear 
the shape of a drum?", and pose the question, "Can one see the hue of the 
Riemann hypothesis?" 

Finally, we turn to a description of the organization of the paper. In section 
1 the problem will be introduced in a form that is conducive to the employ- 
ment of tools which are needed to solve it. Based on three propositions, two 
of which are known results, while the third one still needs to be established, 
the short proof of the major result will be given in that section. In section 2 
material from this author's previous work will be assembled in a fashion that 
facilitates its usage in the present context, and a number of preparatory results 
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will be established. In section 3 the proof for the outstanding proposition will 
be provided. 



Section 1 

The observation that gaps tend to open up readily for small coupling constants 
due to basic perturbations of the degenerate case, that is when the coupling con- 
stant is equal to zero, naturally leads to a search for an argument that allows 
one to show that those gaps can not close as the coupling constant increases. 
The advantage of such an approach goes beyond mere expediency. By estab- 
lishing the persistence of gaps as opposed to proving the existence of gaps for 
specific parameters, one can reach through to the elusive self-dual case "from 
the outside". It is the purpose of the present work to develop such an argu- 
ment. Earlier attempts can be readily traced in the literature. For instance the 
major thrust in [HKS] is along the same line: If one substitutes the conjecture 
C3 in that paper, which at this point remains unproven, by the combination of 
Proposition 2 and Proposition 3 below, then one can simply follow the argu- 
ment provided in "Remarque 3.2.4" in [HKS] to obtain the desired result. In 
conclusion our proof will require two legs to stand on: The first leg establishes 
that gaps with prescribed labels open up for sufficiently small (but by no means 
uniform!) coupling constants. This has been rigorously shown in [HKS]. The 
second leg guarantees that (open) gaps can't close, as the coupling constant 
increases. Our argument to accomplish this has two components. First, we 
need to show that if a gap closes, then it follows that the Lyapunov exponent, 
considered as a function of two real parameters, namely the coupling constant 
and the spectral parameter, has a local maximum in one of the resolvent sets. 
Second, we need to prove that the Lyapunov exponent does not have any critical 
points. It is the proof of the second component that will occupy the main body 
of the paper. 

We turn now to outlining the setting for the proof of that particular component. 
Let a be an irrational number, let u and v be unitary operators satisfying the 
relation 

(1.1) 

uv = e 2m vu, 

and let A a be the (abstract) C*-algebra generated by u and v. Furthermore, 
let t be the unique tracial state on A a . This is a positive linear functional, 
standardized by setting t(<b) = 1, where <b denotes the unit element in A a , such 
that r(ab) = r(ba) holds for all elements in a,beA a . In this setting we define 
for any positive coupling constant (3 the almost Mathieu operator as follows, 

(1.2) 

h{/3) =u + u* +/3(v + v*). 

As usual, the upper right asterisk denotes the adjoint of an operator. The 
integrated density of states can now be identified with the restriction of the 
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functional r to the abelian C*-algebra generated by h(/3). To obtain the inte- 
grated density of states proper, all one has to do is to represent this restricted 
functional by a probability measure on the spectrum of h{(3). In light of the 
comments made above we can now formulate our first proposition which estab- 
lishes the first leg of the argument. 

1.1 Proposition ([HKS]) For any number re[0, l]nZ+aZ there exists a positive 
number /?o with the property that for any /3e(0, /?o], the operator h{[3) has a gap 
in its spectrum with the label r. More specifically there exists a real number 
sp in the resolvent set of h((3), such that the spectral projection p associated 
with the interval (— oo, sp) has the property r(p) = r. 

We turn now to the second leg. First we need to define the Lyapunov exponent 
in a way that is compatible with the present settings. 

L : R x C -> JR, L(/3, z) = r(log \h(f3) - z\). 

While the operator log \h(/3) — z\ = log[(/i(/3) — z)(h(/3) — z)]i is an element of A a 
for ze<C\Sp(h((3)), Sp(h(f3)) denoting the spectrum of h(0), this is not the case 
for zeSp(h(P)). However, for any complex number this operator is contained 
in L 1 (A, t), the space of "integrable operators" associated with A a and r. By 
virtue of the so-called Thouless formula the number L(j3, z) is seen to coincide 
with the usual definition of the Lyapunov exponent for h{0) at z. We can now 
formulate our next statement. 

1.2 Proposition The function L is jointly continuous is both variables. More- 
over, for (3 < 1, L(0, z) = for every zeSp(h((3)). 

While this follows for rotation numbers satisfying a diophantinc condition from 
the present author's earlier work in conjunction with the semicontinuity of the 
spectrum (for badly approximable numbers see [R4], Proposition 2.12, and for 
sufficiently well approximable numbers see [R2], Corollary 2.3), a proof which 
does not rely on any diophantine condition is implicit in [B J] . Indeed, the crucial 
Proposition 9 towards the end of that paper is valid for any pair of sufficiently 
smooth potentials, not just those which differ by the spectral parameter E. The 
remaining argument then carries through essentially without change. 

The third statement, which will be proved in section 2 and section 3, is as 
follows. 

1.3 Proposition For every ze¥L\Sp(h(/3)), 

(T{(h(p)-z)- 1 ),T((h(0)-z)- 1 v))?(O,O). 

The following theorem, whose proof is the main objective in the sequel, has 
many precursors. The most up-to-date partial result appears to be [AJ], where 
the stated claim is that the strong version of the "Ten Martini Problem" holds 
for a set of badly approximable rotation numbers, and for all coupling constants 
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other than and 1. For a survey of contributions that preceded this partial 
result the reader is referred to that paper. 

1.4 Theorem For every re[0, 1] n {% + a%) there exists seM\Sp(h(p)) such that 
the spectral projection p associated with the interval (— oo, s] has the property 
r(p) = r. 

Proof: By duality it suffices to consider the case < /3 < 1 only. Let re[0, 1] (~l 
(7L + a7L), and choose /?o as in Proposition 1.1. For some /3e(0, (3q) let be as 
in Proposition 1.1. Let f2 r be the connected component of R = {(/?, t)elR + x 
1R/ te]R\Sp(h((3))} containing the point (/3, s^). Since (3 can be chosen arbitrarily 
close to zero, it suffices to show that f2 r contains a point whose first coordinate 
is equal to 1. Suppose this were not true. Then fl r C [0,1] x [—4,4], and 
therefore O r is compact. By Proposition 1.2 the function L is continuous on 
f2 r , and it takes the value zero on the boundary dfl r . It follows that L has a 
local maximum at some point (f3 e , s e ) in fl r . Since L is infinitely diffcrcntiable 
in n r , the gradient of L at (/3 e ,s e ) vanishes. Thus 

°— (/3 e , Se ) = r(( Se -M/3e))- 1 ) = 0, 

81 

— s e ) = T((h(P e ) - Se)- 1 ^ + v*)) = 2r((h(fi e ) - s^v) = 0. 
By Proposition 1.3 this is impossible. 

Remarks: 1) While the first of the two partial derivatives occurring in the proof 
is obvious, the second one warrants a few words of explanation. Differentiating 
the first of the above partial derivatives with respect to f3, by invoking the 
formula 

^(z - hi/3))- 1 = (z- h{(3))-\v + v*)(z h{(3))-\ 

and then antidifferentiating the result with respect to z yields the claimed for- 
mula plus a function which depends on (3 only, f{(3) say. In order to show that 
f{[3) is actually zero, one would simply like to let z approach infinity, because 
^ then approaches zero. This is of course impossible, since z is confined to a 
(bounded) gap. Therefore one has to replace the real parameter z by a complex 
one, z = t + ie, for a small positive number e. Repeating the steps just outlined 
in this particular situation, and then letting e approach zero, yields the claimed 
formula for the partial derivative 

2) It is worthwhile mentioning, that at this stage it is already clear that any 
possible critical point for the function L has to be a local maximum. Indeed, 
taking the second partial derivatives 

|^(/3,z) = -t((z-M/3))- 2 ), 
z) = T ((z h((3))-\v + v*)(z HP))- 1 ) = t((z h((3))- 2 (v + „*)), 
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^±( ( 3,z) = -r([(z-hm- 1 (v + v*)?), 

and applying the Cauchy-Schwarz inequality shows that the determinant of the 
Hessian of the function L is strictly positive. Since the diagonal entries of the 
Hessian are negative numbers, the claim follows. 

Section 2 

We now turn to the expansion and refinement of the settings introduced in 
Section 1. In the sequel we assume throughout that a is a fixed irrational 
number, and that (3 is a positive number less than 1. For p, qtlL, we define the 
standardized monomials 

Wpq = e -P1™i u P v 1^ 

and for ze(C\Sp(h(j3)) , 

Cpq(z) = T((h(f3) - Zy^Wpg. 

The standardization ensures that these numbers are real valued whenever z is 
a real number. The double sequence {c pq (z)} solves the following system linear 
equations for s = z, 

(2.1) 

cosiraq(x p+lyq + x p _i A ) + cos irap(x Ptq+ i + x Ptq -i) = sx pq , 

smiraq(x p+lyq - x p ^ q ) - (3 sinirap(x Pyq+1 - x pA -i) = 0, 

for all p, qeTL, except p = q = 0. We shall refer to the system (2.1) with the case 
p = q = exempted by (2.1)*. 

Before we proceed with our objective, we are going to dwell a little on the linear 
system (2.1). First, if one multiplies the second equation by the imaginary 
unit i, and adds the result to the first equation, then one obtains Harper's 
equation on the two dimensional lattice, which is so common in physics, and 
there are no algebraic complications attached to this equation. By contrast, the 
system (2.1) and its truncated version (2.1)*, combines two features which call 
for a specifically designed approach to generate and analyze its solutions (see 
[R3]). On the one hand, the system is largely overdetermined, giving rise to 
redundancies: Asymptotically there are roughly twice as many equations than 
variables. On the other hand, the system is degenerate along the diagonals p = q 
and p = — q: A recursion involving 4x4 matrices to generate the solutions of this 
system collapses as one tries to cross either one of those two axes. The tenet 
underlying the proof of Proposition 1.3 is that these two features, reflecting 
intrinsic properties of the operator h(/3), must also hold the key to its more 
elusive spectral properties. Generally speaking, for numbers s in the spectrum 
of h{fi) the system (2.1) yields uniformly bounded solutions, which are obtained 
by evaluating certain states defined on the C*-algebra A a at the standardized 
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monomials w pq . (In [Rl], where it was shown that the linear dimension of the 
space of uniformly bounded solutions is always either equal to one or to two, 
they were referred to as "eigenstates" .) For s=z in the resolvent set of h(/3), 
the double sequence {c pq (z)} is always exponentially decaying as \p\ — > oo and 
\q\ — > cxo. This is the crucial property we shall exploit in the proof. However, in 
oder to take advantage of it, we first need to "homogenize" the double sequence. 
In other words, we have to find a double sequence solving the homogeneous 
system (2.1), which preserves some measure of that exponential decay, but also 
shares the vanishing conditions that hail from a possible critical point for the 
function L. To this effect we need to delve a bit deeper into system (2.1)*. 

The solutions to this system form a linear space of dimension 6. Up to a scaling 
factor to be determined below, there are four solutions with the property that 
each of those has non- vanishing coefficients only in exactly one of the four sectors 
separated by the lines p — q and p = — q in the two dimensional lattice. The 
components of these four solutions are nothing but the "Fourier coefficients" of 
the resolvent of perturbations of h(/3) which are obtained by multiplying the 
generators u and v with suitable complex numbers of modulus larger than one 
or less than one. Put in technical terms, expanding these resolvents in terms of 
the standardized monomials w pq yields a "non-commutative" multiple Laurent 
series whose coefficients are exactly those solutions. To assign a symbol to 
each of the four solutions, let's say that R p 1 q °\s) vanishes for p ^ 0, R pq 10 \s) 
vanishes for p > 0,R v ° q ' 1 \s) vanishes for q < 0, and R p ° q ^(s) vanishes for 
q > 0. All four of these solutions can be computed by means of a two component 
recursion of the type (3.13) in [R3]. Since we assumed /3 to be less than one, it 
follows that _R( 1,0 )(s) and i?( _1 '°)(s) decay exponentially of order (3 along any 
line in the lattice with slope 1 or — 1, and that R t ^ ) ' 1 \s) as well as R^^ 1 \s) 
grow exponentially of order along the four lines with slope 1 or —1 through 
the points (0,1) or (0, —1) in exactly one direction. Taking the arithmetic mean 
of R^' \z) and R^ 1 '°' l (z), both suitably scaled, yields a solution {d pq (z)} of 
(2.1)* which has the following properties 

(2.2) 

d m {z) = d| p ||g|(z),fora]lp,qreZ;d 1)0 (z) = ^;d pq (z) = Ofor \q\ > \p\ . 

(2.3) 

lim /3~ |p| \d p .k+ P (z)\ < oo, lim \d p ^- P {z)\ < oo, for all ktTL. 

|p|^oc |p|— >oo 

Returning to our objective, it now follows that the double sequence <f> pq (z) = 
c pqi z ) — d pq {z) solves the system (2.1) and has the following additional proper- 
ties: 

(2.4) 

{4>pq( z )} decays two — sided exponentially along any line in the lattice with slope 1 or — 1 
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(2.5) 

If coo (z) = c i(z) = 0, then <j> pq {z) = Ofor \p\, \q\ < 1 

We are going to shelve this for a while, and turn to the discussion of a cer- 
tain subalgebra of the C*-algebra A a . For two elements a, beA a we denote by 
dig* {a, b) the *-algebra generated by a and b. Next wc define two distinguished 
elements. 

U = (3-^11 + f3^v,V = tu_i,i 

which satisfy the relations 
(2.6) 

U¥ = A" 2 ¥U,U*¥ = A 2 ¥U*,whereA = e™ % . 

One way to look at these two elements is that, while mimicking the generators 
for the rotation algebra A a , they also allow for the representation of the element 
h((3) in a form that resembles representing the degenerate element h(0) in terms 
of u and v, 

h(/3) = V + V*. 

The next step is to complete the mimicry, rendering alg*(V, ¥) as much as 
possible a look alike of alg*(u,v). What's missing from the picture is a basic 
symmetry, a conjugate linear involutive automorphism that assigns to one of the 
generators its adjoint, while fixing the other one. Such a symmetry is readily 
available for alg*(u,v), 

(2.7) 

o{u) = u* , cr(v) = v. 

Reversal of the roles of u and v leads to another symmetry, which is equivalent 
to the one just defined. But due to the asymmetric nature of the elements U 
and ¥, only one of them survives the mimicry. To obtain such a symmetry for 
alg*(U, ¥), we first introduce an automorphism of A a that appeared for the 
first time in [R4], 

Pp{u) = vuv(uv + (3) 1 (v*u*+f3) 

pp{v) = v(uv + /3) _1 (t/V + P). 

Notice that, due to the general properties of A a , it is quite easy to define 
automorphisms of A a . All one has to do is to assign unitary elements in A a to 
the two generators which preserve the fundamental commutation relation (1.1). 
Any assignment of this kind extends automatically to an automorphism. Since 
it can be shown that U is a generator for A a , in other words the set of all (non- 
commutative) polynomials in U and U* is norm dense in A a , the automorphism 
P/3 is uniquely determined by the identity, 

(2.9) 

Pfi{u + (3v) = u* + 0v. 
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Thus, composition of a and pp, 

op = a opp, 

yields a conjugate linear automorphism which is uniquely determined by the 
assignments 

(2.10) 

ap(u) = vu*v(u*v + /?) 1 (v*u + f3), 
ap(v) = v(u*v + f3y 1 (v*u + (3). 

The restriction of this symmetry to the algebra alg*(V,V) is exactly what we 
need, 

(2.11) 

ap(V) = V,ap (¥) = ¥*. 

While the first relation is obvious, the second one can be checked through 
straightforward manipulations. Our next objective is to show that this sym- 
metry, when evaluated at the standardized monomials, yields elements whose 
expansion in the standardized monomials have desirable exponential decay prop- 
erties. First we observe that 

oc 

(u*v + (3)- 1 = ^2(-/3) n (v*u) n+ \ 

n=0 

and taking the adjoint on both sides yields of course a similar expansion. This 
shows that ap(w pq ) is a product of v p+q and an clement that has a power series 
expansion in the monomial u>i,_i whose radius of convergence is equal to 
In conclusion, we obtain the following representation, 

(2.12) 

<Jp(w pq ) = V r<£' q) w m , p+q - m , where lim | r ( ^' q) \T^ < (3 . 

In preparation for the proof of Proposition 1.3 in the next section we introduce 
two linear functional (p z and ip z • op which are defined for ze]R\Sp(h(f3)) on 
the algebra alg*(u,v) by the assignments 

Vz(w pq ) = 4> pq {z), 

lf z • (Tp(w pq ) = ^2 r rn q)( Pz(w m , P + q -rn) 
mem 

Notice that, by (2.4) and (2.12) the terms in the sum on the right-hand side of 
the second formula decay exponentially of an order less than or equal to f3. The 
second formula defines essentially the composition of the first functional with 
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the symmetry a p. Since 073 is conjugate linear, we have to conjugate the terms 
in the sum in order to render the resulting functional linear. We are now going 
to show that the two functionals are actually equal. 

(2.13) 

To see this, we need to return briefly to the settings at the beginning of the 
present section. First, since op(h(f3)) = h(f3), we also have op((h(/3) — z)^ 1 ) = 
(h((3) — z) . Notice that, by our assumption, z is a real number. This means 
that all we need to show is, that the functionals $ z and $ z »ap, which are defined 
below, are equal. 

(2.14) 

■&z(w pq ) = R p 1 q ' 0) (z),-& z »ap(w pq ) = r< £' q) ^z(w m , p+q - m ). 

maTL 

and that a similar statement holds for R(~ 1,0 \z). Again, since ap(h((3)) = h{(3), 
and since R^^(z) solves the system (2.1)*, the second functional solves the 
system (2.1)* as well. However, since by the representation in (2.12) and the 
vanishing properties of R^^(z), d z • ap{w pq ) vanishes for all indices located 
below or on the the line p = — q in the two dimensional lattice, the double 
sequence {$ z • ap(w pq )} must be a linear combination of R( 1,0 \z) and R^ 0,1 \z). 
Since _R( 0,1 )(z) grows exponentially of order (3~ 1 along the two lines with slope 
1 or— 1 through the point (0,1), while -d z • op{w pq ) and R p \'°\z) vanish or 
decay exponentially along these two lines, it follows that {i? z • ap(w pq )} is just 
a scalar multiple of i?( 1,0 ^(z). But since op is unital, the two double sequences 
must actually be equal. A similar statement can now be obtained for R^- X ^{z) 
along a similar line of reasoning. This concludes our argument establishing the 
validity of (2.13). 



Section 3 

We turn now to the proof of Proposition 1.3. Henceforth we shall simply write ip 
for the functional ip z introduced in section 2, because we shall assume that z is a 
fixed real number in the resolvent set of h(/3). The idea of the proof is to exploit 
the decay conditions of ip(w pq ) along lines with slope — lin the two dimensional 
lattice, to construct a functional ip on alg*(u,v), with the property that ip can 
be recovered from ip by the identity ip(a\J) — <p(a) for all aealg*(u, v), and then 
showing that such a functional can not exist, in case the function L in section 
lhas a critical point. Of course, there are infinitely many functionals with this 
property. One simply has to implement a separate elementary recursion along 
every single line with slope — 1. The crux is to impose constraints which 
restrict the availability of such functionals severely. More specifically, we shall 
prove the following. 
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3.1 Lemma There exists a linear functional ip on alg*(u,v), having the prop- 
erties, 

(i)Vaealg*(u,v) : ?/>(aU) = <p(a). 



(ii)VpeN ,V<?eZ : ip(JJ p V q ) = V(Up¥"«), ^((U*)^ 9 ) = V((U*)n/-«). 
(iii)Vaea^*(u, v) : V((/i(/3) - z)a) = i>{a{h{(3) - z)U) = 0. 

Proof: First we extend slightly the settings of Section 2. Let 03 be the set of 
elements a in A a which can be written in the form 

(3.1) 

N 

a = k mtn - m w m ,n~m, where lim \k m ,n-m\ |m| < oofor-A^ <n<N. 

' '-^ \m\ >oc 

n=-N mt% 

This set is an involutive subalgebra of A a . Also, (2.12) implies 
(3.2) 

a/j (<8) = 58. 

Furthermore, 

oo 

U" 1 = /H ^(-/3)"(u*w)"ii*e ! B. 

Now (2.3) allows us to extend the definition of the functional tp to elements of 
the form (3.1) as follows, 

JV 



^( a ) — ^ ^ ^ ^ km.n—m^i^m.n—vn)- 



Obviously, this extended functional is also linear. Moreover, by (2.13) 

(3-3) 

y(ap{a)) = </?(a),ae<8, 

and since the double sequence {4> pq (z)} in Section 2 solves the system (2.1), 
(3.4) 

<p{{h{/3) - z)a) = ip{a{h{(3) - z)) = 0, aeQ3. 
We are now going to define the functional ip as follows, 

ip(a) = ip(aV~ 1 ),ae ( S. 

We need to check that ip has the claimed properties. By construction this is 
obvious for (i). Next, if p is a non-negative integer, and q is an arbitrary integer, 
then (2.11) and (3.3) yield, 



ip{V p V q ) = (^(U P ¥ 9 U _1 ) = 1/9(073 (UPVULJ- 1 )) = (^(UPV-sTLJ- 1 ) = ip(VPV-i), 
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which establishes the first identity in (ii). The second identity can be shown in 
the same way. Finally, (3.4) yields, 

i/)((h(0) - z)a) = ip((h(0) - zjar 1 ) = 0, aeQ3, 

and also, 

ip{a{h{/3) - z)U) = f{a{h{/3) - zjUIT 1 ) = <p(a(h(0) - z)) = 0, 
which establishes (iii) as well. < 

Remark: Tracking the significance of the decay condition (2.4) through the 
discussion so far, one observes that this condition is far stronger than what is 
needed to make the arguments work. It would be enough to assume that the 
double sequence in (2.4) does not increase exponentially of an order larger than 
or equal to along any line with slope —1 in the two dimensional lattice. 
All one has to do is to impose a stronger exponential decay condition on the 
elements in the algebra <8. 

In the proof of the following lemma we shall be using nothing but the relation 
(2.6), as well as the following 

(3.5) 

U*U = AV + A~ x ¥* + 7<e, where7 = (3 + (3~ l 

2.2 Lemma If ^ is a linear functional defined on alg*(U, ¥) such that for some 
non-zero real number t the following two conditions hold, 

(i)\/aealg*(U, ¥) : ip((TU + U* - t<e)a) = 0, V>(a(U + U* - te)U) = 0, 

(ii)V(U) = V(U*) = V(¥) = ip{V*) = V(<b) = 0, 
then V(U¥) = 0. 
Proof: Since 

(U + U* - te)U¥« = U 2 ¥ 9 + U*U¥ 9 - WV q , 
¥ 9 (U 2 + U*U - W) = A 4<? U 2 ¥ 9 + U*U¥ 9 - X 2q WV q , 

(i) yields 

^(U*U¥ 9 ) = ^(U¥«) - V(U 2 ¥' 3 ), 
i/j(V*W q ) = t\ 2q ij(W q ) - A 4 V(U 2 ¥ 9 ), 

hence, 
(3.6) 

#(U¥ 9 ) = (l + A 2g )V(U 2 ¥ 9 ),for<7^0. 

Next we derive two elementary identities involving ip(UV) and f/'(U¥ 2 ). Em- 
ploying the second identity in (i) with a = ¥U*, yields 

V>(¥U*(U 2 + U*U - W)) = 0, 
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hence 

V>(V(U*U)U) + V(¥U*(U*U)) - #(¥(U*U)) = o, 
which by virtue of (3.5) yields 

V'(¥(A¥+A- 1 ¥*+7(e)U)+V'(¥U*(A¥+A" 1 ¥*+7(e))-#(¥(A¥+A- 1 ¥*+7(e)) = 0, 
and so by (ii), 

AV>(¥ 2 U) + 7^(¥U) + 7V>(¥U*¥) + 7^(¥U*) - tA?A(V*) = 0, 
which finally implies by (2.6) 

A 5 V>(U¥ 2 ) + 7 A 2 VU¥) + AV>(U*¥ 2 ) + 7A" 2 ^(U*¥) - iA^(¥ 2 ) = 0. 
Combining this with 

V>((U + U*)¥) = #(¥) = 0, 
which follows from the first identity in (i) for a = ¥, and from (ii), we conclude 

A 5 V(U¥ 2 ) + 7 (A 2 - A~ 2 )V>(U¥) + AV>(U*¥ 2 ) - t\^(V 2 ) = 0. 
Combining this with 

V>((U + U*)¥ 2 ) = #(¥ 2 ), 
which is true by the first identity in (i) for a = V 2 , yields 
(3.7) 

A 5 t/>(U¥ 2 ) + 7 (A 2 - A- 2 )V>(U¥) + A(t?/>(V 2 ) - V(U¥ 2 )) - iA^(¥ 2 ) = 0. 
Next, since 

?K(U + U* - te)U¥) = 0, 

which holds by the first identity in (i) for a = U¥, properties (3.5) and (ii) 
yield 

V>(U 2 ¥) + A?/>(¥ 2 ) - #(U¥) = 0. 
Combining this with (3.6) for q = 1 yields, 

—^-2 V(U¥) + AV>(¥ 2 ) - #(U¥) = 0, 
1 + A 

which in turn simplifies to 

^(¥ 2 ) = y-^VtW). 

Combining this with (3.7) yields 

A 5 7/.(U¥ 2 )+ 7 (A 2 -A- 2 )^(U¥)+ T ^V'(U¥)-A^(U¥ 2 )- T ^V'(U¥) - 0, 
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which simplifies to 
(3.8) 

(A 5 - A)V>(UV 2 ) + 7 (A 2 - A" 2 )V(U¥) = 0. 
Next, starting over again, using the first identity in (i) with a = U 2 ¥, 

V>((U + U* -ice)U 2 ¥) =0, 

or equivalently, 

^(U 3 ¥) + V((U*U)U¥) - #(U 2 ¥) = 0, 
yielding together with (3.5), 

i/>(U 3 ¥) + ^((A¥ + A _1 V* + 7<e)U¥) - #(U 2 ¥) = 0, 
which by (ii) implies, 

V>(U 3 ¥) + A 3 V(U¥ 2 ) + 7^(U¥) - #(U 2 ¥) = 0. 
Employing (3.6) to this for q = 1, we obtain, 
(3.9) 

t 2 

V>(U 3 ¥) + A 3 V>(U¥ 2 ) + (7 - —)ip(w) = 0. 

1 + \ z 

On the other hand, the second identity in (i) with a = U¥ yields, 
V>(U¥(U 2 + U*U - tU)) - 0, 

or equivalently 

V>(U¥U 2 ) + V(U¥(U*U)) - #(U¥U) = 0, 

hence by (3.5), 

V>(U¥U 2 ) + V(U¥(A¥ + A _1 ¥* + 753)) - #(U¥U) = 0, 
which by (ii) and (2.6) yields, 

AV(U 3 ¥) + AV>(U¥ 2 ) + 7^(U¥) - U 2 V>(U 2 ¥) = 0, 
and after invoking (3.6) for q = 1 again, 

AV(U 3 V) + A^(U¥ 2 ) + (7 - -^-)^(U¥) = 0. 

1 + A 

Finally, multiplying this by A -4 and subtracting it from (3.9) we obtain, 
(3.10) 

(A 3 - A- 3 )^(U¥ 2 ) + [(1 - A- 4 ) 7 - ^^i 2 ]V(U¥) = 0. 
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Now suppose that 0. Then comparison of (3.8) with (3.10) yields, 

(A 3 - A- 3 )(A 2 - A- 2 ) 7 = (A 5 - A)[(l - A- 4 ) 7 - ^J* 2 ], 
which turns into, 

(A- 6 - A" 4 - A" 2 + 1) 7 = (2 - A 2 - A- 2 )i 2 . 

Since the number on the right-hand side of this equation is real, it follows that 
A~ 6 — A~ 4 — A~ 2 must be real as well or equivalently 

A" 6 - A" 4 - A- 2 =A 6 - A 4 - A 2 

This in turn implies that 

A 12 - A 10 - A 8 + A 4 + A 2 - 1 = 0. 

This conflicts with the fact that, a being irrational, the set {A 2n /neZ} is dense 
in the unit circle. Therefore, ip(JJV) — 0, as claimed. A 

Remark: Performing the kind of manipulations in the proof of Lemma 2.2 for 
more general terms of the form WV q and (U* p )'¥ 9 one can actually show that 
the functional tp "almost" vanishes. If the condition ^(U 2 )=0 is added, then if) 
vanishes completely. 

Proof of Proposition 1.3: First notice that the functional tp can not be 
zero. This is true because the double sequence {c vq {z)} decays exponentially 
as \p\ — ► oo and \q\ — > oo, while this is not true for {d pq (z)}. If {d pq (z)} were 
decaying exponentially as \p\ — > oo and \q\ — > oo, then the two solutions Rifl{z) 
and R-ifi(z) of the system (2.1)* would give rise to two distinct inverses of the 
same element h(/3) — z, which is of course impossible. Now suppose that 

{r((h{0)-z)- 1 ) = r((h{0)-z)- 1 v))=O. 

It follows from (2.5), 

<p(u) = (p(u*) = ip(v) = ip(v*) = 

Since the double sequence {<p(w pq )} solves the system (2.1), this implies that <p 
is non-zero if and only if 

<p(V) ± 0. 

The vanishing conditions for tp translate into several vanishing conditions for ip. 
First, by Lemma 2.1, 

^(U) = <p(<B) = 0. 

Furthermore, 

P~^ip{<R) + (3?ip{u*v) = tp{u*V) = ip{u*) = 0, 
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or, equivalently, 



I3~?ip(v*u) + /3*^(<b) = tp(v*U) = <p{v*) = 0, 
/3~2A^(¥) + / 3^(e) = 0. 



By Lemma 2.1(ii), 

This yields a linear system for ^>(<b) and ip{V), 

/3-iip(<B)+(3i\ip(V) = 0, 

P?iP{<b) +/3-H^(V) = 0. 
Since the determinant of this system, 



(3-2 f32\ 



A(/r 1 - 0) 



is non-zero for (3 ^ 1, it follows that 

ip(e) = tp(¥) = ip(V*) = 0. 

Finally, by Lemma 3.1(iii), since ^(ce) = 0, and V'(U) = <p(<b) — 

0(U*) = ipQU + U* - zee) = 0. 

In conclusion we have shown, that -0 (more precisely its restriction to alg* (U, V)) 
satisfies the conditions (i) and (ii) of Lemma 2.2, letting t = z(3~z. Therefore, 

</j(¥) = ^(¥U) = A 2 ^(U¥) = 0. 

Since we observed at the beginning of the proof that tp(V) cannot be zero, we 
have reached a contradiction. -4 
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